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Abstract A simple, one-equation, galaxy formation model is applied to both the 
halo merger tree derived from a high-resolution dissipationless cosmological simu- 
lation and to a quarter million Monte-Carlo halo merger trees. The galaxy formation 
model involves a sharp entropy barrier against the accretion of gas onto low-mass 
halos, the shock heating of infalling gas far from the central regions of massive ha- 
los, and supernova feedback that drives the gas out of shallow halo potential wells. 
With the first approach, we show that the large majority of galaxies within group- 
and cluster-mass halos, known to be mainly dwarf ellipticals, have acquired the bulk 
of their stellar mass through gas accretion and not via galaxy mergers. With the sec- 
ond approach, we qualitatively reproduce the downsizing trend of greater ages at 
greater masses in stars and predict an upsizing trend of greater ages as one proceeds 
to masses lower than IO'^Mq. We find that the fraction of galaxies with very young 
stellar populations (more than half the stellar mass formed within the last 1.5 Gyr) 
is a function of present-day stellar mass, which peaks at 0.5% at merit =10^'^~^'^Mq, 
roughly corresponding to the masses of blue compact dwarfs. We predict that the 
stellar mass function of galaxies should not show a maximum at nistars > 10^'^ Mq, 
with a power-law stellar mass function with slope -1.6 if the IGM temperature in 
the outskirts of halos before reionization is set by molecular Hydrogen cooling. We 
speculate on the nature of the lowest mass galaxies. 
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1 Introduction 
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The mass growth of galaxies can occur either by accretion of gas that cools to form 
molecular clouds in which stars form or by galaxy mergers. While spiral disks are 
believed to form through the first mode, it is still unclear whether elliptical galaxies 
are built by mergers or not. We use a very simple toy model of galaxy formation 
101 run on top of the dark matter halo merger tree obtained from a high-resolution 
dissipationless cosmological simulation (CS) to understand how dwarf galaxies ac- 
quire their mass. We also use our model to predict the frequency of galaxies such 
as the very metal -poor (l/50th solar metallici ty) gala xy I Zw 18, for which the bulk 
of the stellar mass is younger than < 1 Gyr ( |[l5l;]27ll . see also lU using HST color- 
magnitude diagrams) or 500 Myr ( ||22[ 11411 from photometric studies). 



2 Galaxy formation model 

Galaxies form in DM halos, and our toy model gives the mass in stars and cold gas, 
m, as a function of halo mass M and epoch z, taking into account the fact that for 
stars to form one needs: 1) gas accretion, which is fully quenched for low-mass halos 
1 261 121; 2) in cold form, which becomes inefficient in high-mass halos iH IitII : and 



3) to retain the interstellar gas against supernova (SN) winds 1@1: 

fjkjf \ circ reion 7b^" 
«stars(M,z) = — — ^ — — — , (1) 

where fo-Q^^jQ^ ^ 0.17 is the cosmic baryon fraction, Ureion is the minimum halo 
circular velocity for gas accretion (which rises abruptly after reionization; nJstars = 
for ^circ < t'reion)- i^SN is a characteristic velocity for SN feedback, Mohock represents 
the transition from pure cold to mainly hot accretion, and g is the fraction of retained 
baryons in the form of cold gas. 

Figure[T|describes the efficiency of galaxy formation, 'Wstars/(/b^)- at ^=0 and 5, 
with the parameters tuned to match the observed (z^O.l) galaxy stellar mass func- 
tion (MF) of jstl. Galaxy formation occurs in a fairly narrow range of halo masses, 
that varies with redshift. This is further illustrated in Figure |2] which shows the 
time evolution (from right to left) of halos ending with different circular velocity at 
Z-Q- Interestingly, halos evolve, on average, with nearly constant circular velocity 
(however individual halos have constant M, hence decreasing v^irc in their quiescent 
mode, so this decrease is compensated by the increases caused by mergers). Also, 
the temperatures affected by thermal instability (TI) are fairly narrow for the lowest 
mass galaxies, and TI can therefore be neglected to first order in our analysis. 



We first apply equation 
halos (AHF algorithm of 



to the merger tree obtained from the halos and sub- 
1) of a high resolution dark matter CS. When a halo 
enters a more massive one, it becomes a subhalo and its galaxy becomes a satellite. 
The subhalo orbit is followed until either 1) dynamical friction (DF) causes it to 
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Fig. 1 Illustration 0] of the toy model of galaxy formation (eq. (T) with g = 0, t)ieion=40kms~', 
i;sN = 120kms^', Mshock=8 x 10"/!~'Mo and 40% of the stellar mass stripped at every passage 
through a parent halo) at z=0 (upper halo mass scale) and z=5 (lower halo mass scale). The varia- 
tions with redshift come from the redshift modulation of the dependence of Ucirc on M. The dashed 
cun'e shows the effect of our iterative correction for g (§[3). The shaded region comes from abun- 
dance matching flsll . The disagreement, at the high end, of our model with the abundance matching 
prediction suggests that galaxy mergers are important at high halo masses (see 



fall to the halo center and necessarily see its galaxy merge with the central galaxy 
of the parent halo, or 2) it is tidally stripped and heated by the global halo poten- 
tial to the point that there are insufficient particles to follow it. In this latter case 
(^subhaio < 1.5xlO^/i 'Mq, Corresponding to Vcuc — 17kms-i atz = 0), we assume 
that the satellite galaxy merges with the central one after a delay set by DF, for which 
we adopt the timescale carefully calibrated by lll6ll with hydrodynamical CSs. 

The present-day relation between stellar and halo mass obtained when our model 
(with g - 0, which is adequate for massive galaxies) is applied to the CS fl\ is 
impressively close to the prediction derived by abundance matching Ll3i1 . Moreover, 
the variation of the stellar mass functions in different bins of halo mass matches very 
well the measurements from the SDSS i28l|. 

Figure |3] shows that while mergers dominate the growth of the massive galaxies 
(as expected from the toy model, since gas accretion is quenched at high masses), 
their importance drops sharply when one moves to stellar masses below 10''/riMs 
(the mass resolution is ni^^ ^ 10"^'^/!~'Mo, where the median fraction of stellar 
mass acquired by mergers no longer decreases with mstars faster than in our reference 
model where baryons trace the dark matter [red dotted line]). This dominance of gas 
accretion at low mass is also true for the satellites of clusters (dashed line). Since 
observations indicate that most satellites of clusters are dwarf ellipticals (dEs), we 
conclude that cluster dEs are most often not built by mergers. One must resort to 
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Fig. 2 Halo circular velocity versus redshift. Logarithmic mean evolution {solid curves); typical 1 
and 2o" fluctuations in ACDM scenario (dashed curves). The thick broken solid line indicates our 
adopted history of the temperature of the IGM (just outside the viiial radius of galaxies), while the 
cyan hashed region around it shows the uncertainty on this temperature. The equivalent circular 
velocity is the minimum for gas accretion onto halos. The thick rising magenta curve is the limiting 
circular velocity for shock heating of infalling gas near the virial radius. The gray horizontal bands 
indicate temperatures where solar and subsolar (narrow hashed) and also 1% solar (wide hashed) 
metallicity gas is thermally unstable. 



Other mechanisms (not included in our toy model) that transform dwarf irregulars 
into dEs (e.g. harassment ll20ll or ram pressure stripping 



3 When do dwarf galaxies form their stars? 



We have used the halo merger tree code of [21] to statistically study the star forma- 
tion histories (SFHs) of dwarf galaxies. We consider 24 final halo masses geometri- 
cally spaced between 10^ and IO^^'^^/i'^Mq, and run each halo merger tree 10000 
times. We run the toy model (eq. HI with = 0) on the branches of the halo merger 
tree (moving forward in time) to follow the evolution of stellar mass. To compare 
our predictions to SDSS measurements (at z - 0), we correct our z-Q stellar masses 
for the non-inclusion of gas in equation (HJ, assuming that the gas-to-star fraction is 
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Fig. 3 Median fraction of z=0 stellar mass (for /imstars > 10'"Mo one can neglect the gas) acquired 
by mergers, for all galaxies {solid green) and ftMhaio > IO'-'Mq cluster satellites (dashed black). 
The red dotted line shows the baryons trace the dark matter model. The enm bars are uncertainties 
on the median from 100 bootstraps. 



— = = dex {a + b\og mstars ) , (2) 

S \-g 

where we derive a - 4.5 and b - -0.5 by fitting a straight Hne through Fig. 11 of 101 ■ 
This correction involves solving an implicit equation, which we perform iteratively. 

The predicted MFs match fairly well the observed MPs. The left panel of Fig.|4] 
suggests that reionization must occur late (z=6). The middle panel hints that, be- 




Fig. 4 Etfects of thermal history of IGM on galaxy stellar mass function. Left: Effect of reioniza- 
tion epoch. Middle: Effect of pre -reionization IGM temperature. Right: Effect of post-reionization 
IGM temperature. In all plots, the symbols represent the observed SDSS stellar mass functions, 
measured by (sj] (purple asterisks), U (brown circles, strongly incomplete below logms,ajs = 7.4), 
and lisll (gray triangles). 
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fore reionization, the temperature of the IGM must be set by molecular cool- 
ing (ypre-reion - 2kms~'). The right panel indicates a good match between pre- 
dicted and observed stellar MFs when the IGM after reionization is not too cool 

(Vst-reion>30kms-'). 

According to the left panel of Figure |5] as one proceeds from the highest final 
stellar masses to lower ones, the median stellar age first diminishes, qualitatively 
reproducing the downsizing of star formation. However this downsizing stops at 
wJstare ~ 10"^Mo and as one proceeds to even lower masses, one notices an upsiz- 
ing of stellar ages, first weak, becoming strong at mstars < IO^'^Mq. In our model, 
the smallest galaxies form most of their stars before reionization. Classical dwarf 
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Fig. 5 Left: Median (mass-weighted) star fonnation redsiiift vs z=0 mass. Tlie points ( 1 in 5 plotted 
for clarity) are individual galaxies, while the red symbols are medians (error bars extend from 
16th to 84th percentiles). The stripes are artifacts of our discrete set of final halo masses. Right: 
Fraction of young galaxies (half the mass fornied within last 1.5 Gyr) versus z=0 stellar mass (with 
numbers of expected SDSS galaxies). In both plots, we adopt t)pre-reion=2kms~' (H2 cooling) 
before reionization (z=6), Upost-reion^SOkms"' forz<6, usN=120kms~', and /iMsi,ock=8xlO" Mq. 



spheroidals, with mstars = IO^Mq should have moderately old stellar populations, 
with a large scatter (caused by the nearly constant average evolution ofhalo circular 
velocity), while ultra-faint dwarfs should be extremely old (see also ||23J). Note that 
our model predicts a discontinuity in median age versus mass at x lO^M©, which 
appears to be in conflict with the continuity of the metallicity-mass relation lfl9ll . 

Young galaxies must cross the low-mass (entropy) barrier (thick broken line in 
Fig.© only very late. The right panel of Figure |5]shows that the frequency of galax- 
ies with the bulk of their mass in stars acquired within the last 1 .5 Gyr is maximal 
at ^ 0.5% at mstars = lO^'^~^'^M0. The presence of a young galaxy such as I Zw 18, 
whose stellar mass is of order IO^^Mq 12711 . is consistent with our model; we pre- 
dict 0.02% of young galaxies at \ogm - 7 + 0.5, and with 941 in this mass range in 
the SDSS, we expect 0.2 young galaxies in SDSS at this mass, hence (from Poisson 
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Statistics) there is a probability of 18% of detecting at least one galaxy as young as 
I Zw 18 in this mass bin. 

Our model also predicts several hundred young galaxies in the SDSS, mostly at 
'"stars = IO'^Mq (fraction of 0.2%). In comparison, a first analysis of luminosity- 
weighted ages of high surface brightness SDSS galaxies by |ld] leads to ^ 100 
times higher fractions of young galaxies: 16% with ages younger than 1 Gyr at 
m - IO'^Mq. a second analysis by the same team, this time of mass-weighted 



stellar ages for a large fraction of SDSS galaxies 111 111 , reveals much lower frac- 
tions of young galaxies: 1.4% at log;7Zstars = 10 + 0.5, but still as much as 15% at 
log 'Wstars — 9 + 0.5, Still respectively 7 and 30 times higher than our predictions. 
Also, the fraction of ages (summing over all masses) less than 4 Gyr derived from a 
much more-refined semi-analytical model jKl run on the Millennium cosmological 
dark matter simulation jzsl] is negligible in comparison with |ll]'s prediction of 
3%. This discrepancy in the predicted fractions of young stellar populations among 
galaxies of intermediate mass is thus an open issue worth further exploration. 



4 The lowest mass galaxies 

The middle panel of Figure|4]indicates that there is no peak in the best fitting galaxy 
MF, if the IGM temperature before reionization is set by molecular Hydrogen cool- 
ing (t;pre-reion=2kms~') with a low-end slope of -1.59. This is in excellent agree- 
ment with the slope of -1.58 found by If, instead, the IGM temperature before 
reionization is set by atomic Hydrogen cooling (ypie-reion=17kms"'), the mass func- 
tion peaks at ?Mstars = 

10^-5 /i-i Ms . This maximum is probably not caused by of our 
mass resolution, since no such peak is seen when the pre-reionization IGM temper- 
ature is set by H2 cooling. 

The importance of the low-end tail of the galaxy MF raises the question of the 
nature of very low mass objects (wstars < IO^Mq). Two classes of objects come to 
mind: Globular Clusters (GCs) and High Velocity Clouds (HVCs). However, in our 
model, these objects must be (or have been) associated with DM halos. While Galac- 



tic HVCs do appear to require DM |6i|. Galactic GCs don't (e.g. 12411 '). perhaps be- 
cause they are closer and more tidally stripped. 
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